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Abstract. Although inflationary models generically predict a flat spectrum of gravitational 
waves, we point out a general process that produces a sharply peaked spectrum of gravi- 
tational radiation. This process is generic for inflationary models with a complex inflaton 
field which couples to fermions. In particular, for chaotic models these may be the most 
extreme gravitational waves in the Universe with a very large abundance Qcwh 2 ~ 10~ 9 and 
ultra-high frequency, 10 10 Hz. Although not amenable to space based interferometers, the 
signal from this model may be detectable by future table top experiments. 
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1 Introduction 

In addition to providing solutions to the horizon and flatness problems, inflation is also 
responsible for the generation of density perturbations and reheating the universe. The 
excitation of both scalar and tensor perturbations on very large scales with a nearly flat 
spectrum is a generic consequence of inflation. For the highest scale model of inflation, 
p ~ 10 64 (GeV) 4 , it is possible to obtain a large tensor to scalar ratio, r ~ 0.1 — 0.001, which 
can potentially be observable with the PLANCK satellite [1]. 

Lack of any positive detection of stochastic gravitational waves (GWs) suggests that 
the scale of inflation is lower than the upper limit 10 64 (GeV) 4 . Chaotic inflation with 
super-Planckian VEVs is perhaps one of the best known examples of high scale which is 
compatible with the current CMB data [2]. In this scenario, one has a simple potential given 
by V = m 2 (j) 2 where m ~ 10 13 GeV is the mass of the inflaton field <fi [3, 4]. Inflation occurs 
at super-Planckian VEVs, ~ 5 — 10Mp (where Mp = 2.4 x 10 18 GeV), and ends near 
(f> ~ Afp, and can potentially lead to long wavelength stochastic gravitational waves [1]. 

However, there is another interesting source of generating GWs connected intimately 
with the end of inflation. The end of inflation gives rise to an epoch of reheating which 
can be responsible for generating GW, albeit at short sub-Hubble scales, with a spectrum 
which would be peaked as compared to a flat spectrum. The process of reheating is es- 
sential for generating thermal entropy and the creation of the observed degrees of freedom. 
After inflation, it is typically assumed that the inflaton oscillates coherently and eventu- 
ally decays through ordinary perturbative processes (2 or 3 body decays of the inflaton) or 
non-perturbatively (for a review see [5]). However a novel possibility has been suggested in 
which inflaton decay could be delayed if the inflaton forms non-top ological solitons of its 
own matter: inflatonic Q-balls [6, 7], which subsequently decay via surface evaporation [8]. 
Such fragmentation leads generically to a very inhomogeneous and anisotropic distribution 
of matter. If the matter undergoing fragmentation once dominated the energy density of 
the Universe, then the present-day abundance can be very large. This idea has been ex- 
plored in the context of super symmetric flat direction condensates originally in [9, 10] and 
more recently by [11, 12]. 1 The cosmological implications of super symmetric Q-balls has 

1 The Q-balls are formed from the fragmentation of the inflaton condensate. The fragmentation of any 
scalar condensate was first described in the context of gauge mediated supersymmetry breaking [13, 14] and 
subsequently in gravity-mediated models [15, 16]. In these supersymmetric models the scalars are made up of 
combinations of squarks and sleptons. The potential obtains quantum corrections by virtue of gauge/gaugino 



- 1 - 



been extensively studied [9, 15-18, 20, 23-52]. It is this fragmentation process, leading to 
inhomogenous and anisotropic motion of matter in the scalar condensate that sources the 
GWs. 

In this paper we point out that the fragmentation of an inflaton condensate would 
naturally lead to a very large amplitude GW with a high peak frequency of order ~ 10 10 Hz 
with a narrow spectrum which falls as k~ l on sub-Hubble scales. Such GWs may be detectable 
in foreseeable future table-top experiments [53-57] . Indeed the authors of [57] have obtained 
an impressive direct bound Qcwh 2 < 6 x 10 25 in a narrow bandwidth around 100 MHz. 
We should note however that the peaked GW spectrum predicted in this paper coming from 
inflaton fragmentation can easily occur at lower frequencies for low-scale models of inflation 
and reheating. The only requirements for inflaton fragmentation are that the inflaton be 
complex (carrying a conserved U{\) charge) and that the potential around the era of reheating 
be less steep than <f> 2 . Here we focus on chaotic inflation and leave a detailed survey of other 
models for future work. 

We note one further possible signature of high frequency GWs obtained from cosmic 
microwave background (CMB) measurements. If the strength of primordial magnetic fields 
is sufficiently large, it may be possible to obtain indirect bounds on such GWs from CMB 
data [58] . The observation of such GWs would provide an unprecedented view of the earliest 
moments of the Universe's history and energies well above those likely to be probed at particle 
colliders. 

We begin in Section II with a general discussion of condensate fragmentation and Q- 
ball formation. Although this discussion is general, we specialize to the case of running-mass 
chaotic inflation for the remainder of the paper for ease of calculation. In Section III we 
analyze the rate of Q-ball decay which is an essential parameter governing the process of 
reheating. In Section IV we derive properties of the GW produced from the fragmentation. 
These properties include the peak amplitude and frequency as well as a schematic form of 
the complete GW spectrum. 

2 Instabilities in the condensate 

In order to illustrate this, let us consider a simple toy model of running mass inflaton with a 
negative Logarithmic correction, see for discussion [1]. This may happen if the inflaton decay 
occurs through Yukawa couplings to fermions, i.e. Ci n t ~ hcpipip (where tp is the fermion, 
and h is the Yukawa coupling) 2 . Otherwise the inflaton can have gauge couplings in a 
supersymmetric theory, where the gaugino loops dominate to generate a one-loop correction 

and fermion/sfermion couplings such that the minimum of the potential is flattened and becomes shallower 
than quadratic, for a review see [17, 18]. The scalar condensate oscillating in such a potential gives rise to a 
negative pressure on average [15] , which triggers the instability in the condensate and thus the fragmentation. 
The process of fragmentation has been studied analytically and numerically [19-22]. 

2 Inflaton coupling to another scalar will not lead to a complete decay of the inflaton. However, such 
coupling can trigger a phase of preheating via parametric resonance [59] or a tachyonic instability [60] which 
may produce additional GWs. Here we will consider a minimal scenario for the inflaton to decay to fermions 
and study the consequences. Further note that bounds on the Yukawa coupling to light fermions in running 
mass inflaton models are stringent [61], but weak for heavy fermions. The couplings of the inflaton to these 
fermions need not directly reheat the Standard Model (SM) degrees of freedom as in the case of more realistic 
models of inflation, see [62-64]. One might still require these fermions to eventually decay into the SM degrees 
of freedom which itself is a challenging problem. 
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to the inflaton potential: 



V^)=m 2 \ ( f ) \ 2 i^l+Klog^j +■■■ (2.1) 

where K can be positive or negative. Since in this paper our aim is to set a limit on the 
largest amplitude and frequency GW which can be produced after inflation, we would mainly 
concentrate on the first term during inflaton oscillations. 

The potential Eq. (2.1) is consistent with the COBE normalization of the density per- 
turbations if one takes 5p/p ~ m/Mp ~ 10~ 5 . Implying that the inflaton mass must be 
about m ~ 10 13 GeV as noted above. Note also that quadratic chatoic inflation predicts a 
spectral index n s = 0.97, in close agreement with the current WMAP data [2, 61]. 

Typically, for matter couplings the value of K is negative. As a result when the inflaton 
oscillates coherently, it sees a shallower potential: V ~ m 2 cf) 2 ((f)/Mp) 2K oc (j) 2+2K , which 
results in a negative pressure on average for an oscillating inflaton, ie. p ~ — (\K\/2)p [15]. 
It is the presence of a negative pressure for certain sub-Huble modes that lead to their 
exponential growth and ultimate fragmentation of the condensate. 

Although the process of fragmentation is non-perturbative [10], the largest instability 
mode can be tracked analytically in the linear perturbation analysis of the fluctuations of 
the inflaton condensate. Let us decompose the condensate into a spatially homogeneous part 
and a spacetime dependent fluctuation: <fi — > 4>{t) + 5(j)(x,t) and 9 — > 9(t) + 56(x,t). Then 
the classical equations of motion in a Fridemann-Robertson- Walker (FRW) universe for the 
inhomogeneous part are: 

59 + 3H59 - ^-59 + ^$59 + ^<ty - ^6</> = 0, (2.2) 
a z (f> 9 (p 

V 2 

6<j) + m54> - -=-S<f> - 294>59 + 54> V"{6) - 9 2 5cp = 0, (2.3) 



where we have used V'{4> + 5(f>) = V(4>) + 5(pV"((p), which is accurate in the linear regime: 
5<fi < 4>. 

5<p = 5<p e a{t)+ikx , 59 = 59 e a{t)+ikx . (2.4) 

A solution that has real and positive a{t) is a signal that the condensate is unstable. We can 
apply the ansatz (2.4) to equations (2.2) and (2.3) to find a dispersion for a{t) 

d 4 + ( 2^ + V" + 39 2 ) « 2 + + V" ~ ^) ^ = °. ( 2 - 5 ) 



where we have assumed a a 2 . The positivity of a(t) requires that the last term in Eq. (2.5) 
be negative. Thus we find the instability band to be 



k 2 / ■ \ k 2 

< !L < (9 2 -V") =^f£. (2.6) 



The best-amplified mode h^ est often lies close to the middle of the instability band [19-22]. 
We can maximize a(k) with respect to k to find 

kLt = (V + 79 2 ) = k 2 max Q - «™ , ) (2.7) 
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where in the second line we used the definition of k"^ ax to eliminate V" . 

Using the fact that at an extremum a(k) satisfies a 2 (k) = k\ ax j1 — k 2 , we find the rate 
of growth for the best-amplified mode as 

a(hest) 2 = ^ (2.8) 

V best J ^ V I 

In the limit that the Hubble expansion is negligible one can use the classical equation of 
motion for (j> to show 

YJS. ~ m 2 [i + K + Klog(6 2 /2M 2 )] (2.9) 
9 

When this approximation is valid, one can write all the relevant expressions of the fragmenta- 
tion kbest, k max , and d(kb est ) in terms of V" and V . We list for completeness the parameters 
characterizing the fragmentation: 



u2 = l_ yii k 2 = r 2 I 1 ^max 

^max i v > "'best ""max 



2 16 V 



a{k best ) =^Q~yi 1 • (2-10) 

These expressions emphasize the fact that when Hubble expansion is negligible and nonzero 
angular motion exists, the nature of the fragmentation is determined entirely by the potential 
and its derivatives. Note that when the second derivative of the potential is small V" V/(f>, 

klax « V'/ 9 . 

Let us note that the era of Q-ball formation is bounded by causality. The best-amplified 
mode as given by Eq. (11) for the running- mass potential is k? t ~ m 2 \K\ (1 — |X|/4). For 
this mode to be within the causally connected horizon, the size of the fragmented region 
cannot exceed that of the Hubble length at the time of fragmentation: 

r* 1 < TT- 1/2 m\K\ 1/2 ~2x 10 12 GeV (2.11) 

for in ~ 10 13 GeV and \K\ ~ 0.1. The actual time scale for the formation of lumps depends 
on the growth rate of forming Q-balls, which is roughly given by: 

tJ^^Klm, (2.12) 

where 7 ~ 0.1 — 0.001 estimated from the numerical simulations [6, 10]. The estimated charge 
in each inflatonic Q-ball is roughly given by [6] 

Q = —^—nQ « A 7 2 \K\ 2 rimM 2 

«7x 10 6 . (2.13) 

with 7 ~ 10 -2 . Since the mass of gravity-mediated- like Q-balls scales as M(Q) ~ mQ, these 
are generically unstable. As we will see in the next section, there are a number of processes 
that contribute to total rate of Q-ball decay. 
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3 Decay of Q-balls 



The question of the rate of decay of these Q-ball is crucial to their impact on the GW signal. 
Moreover if the Q-balls do not decay sufficiently fast then they can come to dominate the 
energy density of the Universe. The effect of a period of Q-ball domination on GW signals 
has been studied in the Affleck-Dine case by [11]. Here however we study effects that preclude 
such a period. 

Although we assume that the renormalizable scalar potential is U(l) invariant, there will 
in general be non-renormalizable operators indued by gravitational correction which violate 
this U(l) symmetry. Of course there may be additional sources at low-energy physics with 
U{\) violation, conservatively we can simply take that the scale of such U(l) violation is 
M ~ Mp. Thus we include terms, such as 

(S n + (h* n ) 

VnrM = A ^ M +! 4 \ (3.1) 

Such terms lead to departure from purely radial motion of the complex field <3? = 4>e lS in the 
internal U(l) space, giving the condensate a nonzero £7(1) charge since uq = 0\4>\ 2 [65]. The 
U{1) violating interaction would break the condensate into smaller Q-balls [10]. This process 
can proceed perturbatively [10] as well as non-perturbatively [47]. The lowest order operator 
n = 5 would yield a decay rate, given by: ~ (1 / Q)(dQ J dt) ~ an^ ~ ^ 2 \K\ 2 (fa/Mp) 2 m, 
where a ~ (1/Mp) 2 and ~ ^iK^mcj) 2 . For the largest amplitude of the oscillations, 
fa < Mp, the fastest decay rate of the Q-balls is given by: < ^y 2 \K\ 2 m. This time scale 
is fast but cannot prevent the initial inflaton fragmentation, which leads to the formation of 
Q-balls, i.e. tj 1 2> Tj. 

Furthermore, for a rotating inflaton condensate, there also exits the traditional pertur- 
bative decay channel [65, 66] with a relatively slower rate T ~ h 2 m/8ir. All of these would 
contribute to: lumps of inflatonic Q-balls and relativistic matter. The latter might not have 
come into full thermal equilibrium at initial stages, and would have an initial momentum 
~ m. This ambient plasma while interacting with Q-balls also would trigger rapid evapora- 
tion and diffusion of the Q-quanta from the surface area of the Q-ball. The diffusion rate is 
typically given by [22]: 

I\k// = -4vraT (3.2) 

where a ~ 4 — 6 is a diffusion coefficient [36]. The diffusion rate turns out to be the fastest 
way to get rid of the Q-balls from the ambient plasma. For temperatures below the initial 
momentum of the relativistic fermions in the ambient plasma, i.e. T < m ~ 10 13 GeV and 
f° r Q £ 3 x 10 9 the diffusion rate leads to a complete evaporation of the Q-balls within 
one oscillation period^ m. The decay of the Q-balls leads to the domination of relativistic 
fermions with an initial momentum of the order of the inflaton mass, which could be regarded 
as the instantaneous temperature of the relativistic plasma, Ti ns t ~ m. 

4 The spectrum of GW 

The time of maximal GW production and Q-ball formation both occur when the fluctuations 
go non-linear 5(fi/4> ~ 1. Note that Q-balls form when Hf ~ \K\^l 2 m and oscillations begin 
around Hi ~ m. We can therefore estimate the decay in the amplitude of the oscillations 
as (fro ~ 4>i(H f I 'Hi), where fa ~ Mp, denotes the inflaton field value at the end of chaotic 
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inflation. This yields (j>o ~ \K\Mp, which since 5(j) ~ </>0) can be used to estimate the GW 
energy density at the time of fragmentation as 

\K\ 2 m 2 ^ 

PGW = TTo ' i 4 ' 1 ) 



M 2 



such that 



n OT{t/) ^ K Wffiy. (4 . 2) 



By the present era this energy density has been substantially redshifted 

2 



«WtoO-^-«W („)(*) (g) 

'300\ 1/3 o / x 



1/3 

io- 5 ?r 2 ( ncwCv). ( 4 - 3 ) 

Thus the present-day energy density fraction of GW from the fragmentation of the inflatonic 
condensate can be as large as Qcwh 2 (to) « 10 -9 . 

The peak frequency of the gravitational wave spectrum will be largest when the con- 
densate lumps are unstable and decay immediately on cosmological time scales. In this case 
the peak frequency at the present epoch is related to the frequency at production time via 
the assumption of adiabatic evolution, a(T) oc g~ l ^T~ l . Numerical and analytic calcula- 
tions [6, 13, 19-22] of the fragmentation process indicate that frequency of the fastest-growing 
mode is /* = A^ est /27r where /~ is the time of fragmentation [10]. Then the frequency 
/o = /* (a*/ao) is 



/o — /* 



9s,o \ 1/3 / Jo 



O'inst J \9s,R ) \Tifist 



assuming a* ~ ai ns t, where all parameters with subscripts "inst" and "0" are evaluated at the 
epoch when the universe became radiation dominated with an instantaneous temperature, 
Ti ns t ~ 77i, and the present epoch respectively. This is by far the highest frequency GW 
which can be generated by any mechanism without affecting the scalar perturbations of the 
cosmic microave background radiation. 

The full spectrum can be estimated using existing methods for GW from first order 
phase transitions, see for instance [67, 68]. In this formalism one begins with the standard 
linearization of the Einstein equation 



32 

Uhij = —nGa 2 pU l:j (4.5) 



where p is the energy density in the fragmenting condensate and Iljj is the anisotropic stress- 
energy tensor. All the relevant information of the source is encoded in the power spectrum 
P s (k,t,t') which can be obtained from the anisotropic tensor via the 2-point correlation 
function 

(n^MlWk'/)) = (2n) 3 5(k-k')P s (k,t,t'). (4.6) 
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Figure 1. Here we plot the GW spectrum from inflaton fragmentation under the assumption that 
the fragmentation is an incoherent process, i.e. the anisotropic stresses are uncorrelatcd at different 
times. 



Using the definition of gravitational wave energy density pcwi*) = (^«j( x )^ij( x ))/(87rGa 2 ) 
it can be shown that the present-day spectrum is given by 



GW 



4Q 



rad,0 



d log k 3-7T 2 



rad 



Htk 6 Re[P s (k,k,k,)], 



(4.7) 



where we have assumed that the GW production occurs during radiation domination. Here 
we make the simplifying assumption that the power spectrum is separable at equal times 
P s (k,t,t) = \F(k)\ 2 \g(t)\ 2 . The time dependence of the fragmentation can be taken to turn 
on abruptly and last for a duration T" 1 such that g(t) is unity during fragmentation and 
zero otherwise. The spatial dependence |-F(A;)| 2 is taken to have the simple form 



\F(k)f 



R 3 



1 + (kRY 



(4.8) 



where R = 2i:k^ est is the characteristic scale of the fragmentation and the rate of frag- 
mentation is approximately T = 7 _1 ff^. Moreover the fragmentation is a temporally random 
process. Thus in our case the GW source is totally incoherent in the sense that the anisotropic 
stress at different times is uncorrelated. Eq. (4.7) is plotted in Fig. 1 under the assumption 
the fragmentation is temporally abrupt, short, and incoherent. 

Before we conclude there are couple of points to note. The abundance of GW in the 
aftermath of inflation is the largest amongst any of the other competing scenarios. Scalar 
preheating can also create large amplitude GW with an abundance Qcwh 2 ~ 10~ 10 , but 
with a peak frequency / ~ 10 7 Hz [59], as compared to the fragmentation of the inflaton 
with Q,G\yh 2 ~ 10~ 9 at peak frequencies around 10 10 Hz. In order to exceed the above limit 
one would require the universe to be filled with large energy densities and larger Hubble 
expansion rate. 

Note that in order to produce the large scale structures and the amplitude of the seed 
perturbations for the temperature anisotropy of order 10 -5 , the inflaton energy density can 
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not be arbitrarily large, as it is bounded by the null observations of the stochastic GW 
which are inevitable consequence of inflation. However if future observations point towards 
peak frequencies in GWs exceeding 10 10 Hz, then this would suggest a sub-Hubble origin for 
GWs at higher Hubble expansion rate and at higher energies then currently envisaged scale 
of inflation which is bounded by 10 64 (GeV) 4 . This would perhaps be the cleanest way to 
falsify inflation as a source for seeding perturbations for the structure formation. 

5 Conclusions 

We have shown a simple and generic process that can produce extreme gravitational radiation 
with a sharply peaked spectrum, in contrast with what one ordinarily expects to produce 
during inflation. The process of reheating in the early Universe is crucial for the production 
of entropy and the observed degrees of freedom. Reheating is however a very difficult era 
to constrain by observation. As we have shown above the if reheating involved an era of 
inflaton fragmentation, then it can have a critical impact on the prospects for GW detection. 
In particular, we have shown that a very early stage of inflaton fragmentation can produce 
extreme GWs with appreciable abundance in the multi-Gigahertz frequency range. Such 
large frequencies make these GWs easily distinguishable from models of Affleck-Dine frag- 
mentation, and the properties of the spectrum shown above make it distinct from what one 
would ordinarily predict from inflation. Although these GW are not amenable to observation 
via space-based interferometers it is conceivable that future tabletop experiments [53, 54] 
may be able to reach the required sensitivity in order to confirm or exclude an era of inflaton 
fragmentation in the early Universe. Moreover, it may be possible to observe signatures from 
such GWs in the cosmic microwave background [58] . Such an observation or refutation would 
give us profound insight into the nature of inflation and the era of reheating. 
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